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ABSTRACT 

We present the results of extensive multi- frequency monitoring of the radio galaxy 3C 111 between 
2004 and 2010 at X-ray (2.4-10 keV), optical {R band), and radio (14.5, 37, and 230 GHz) wave 
bands, as well as multi-epoch imaging with the Very Long Baseline Array (VLBA) at 43 GHz. 
Over the six years of observation, significant dips in the X-ray light curve are followed by ejections 
of bright superluminal knots in the VLBA images. This shows a clear connection between the 
radiative state near the black hole, where the X-rays are produced, and events in the jet. The X-ray 
continuum flux and Fe line intensity are strongly correlated, with a time lag shorter than 90 days 
and consistent with zero. This implies that the Fe line is generated within 90 light-days of the source 
of the X-ray continuum. The power spectral density function of X-ray variations contains a break, 
with steeper slope at shorter timescales. The break timescale of 13lg^ days is commensurate with 
scaling according to the mass of the central black hole based on observations of Seyfert galaxies and 
black hole X-ray binaries (BHXRBs). The data are consistent with the standard paradigm, in which 
the X-rays are predominantly produced by inverse Compton scattering of thermal optical/UV seed 
photons from the accretion disk by a distribution of hot electrons — the corona — situated near the 
disk. Most of the optical emission is generated in the accretion disk due to reprocessing of the X-ray 
emission. The relationships that we have uncovered between the accretion disk and the jet in 3C 111, 
as well as in the FR I radio galaxy 3C 120 in a previous paper, support the paradigm that active 
galactic nuclei and Galactic BHXRBs are fundamentally similar, with characteristic time and size 
scales proportional to the mass of the central black hole. 

Subject headings: galaxies: active — X-rays: galaxies — X-rays: binaries — Galaxies: individual (3C 
111) — quasars: general — Black hole physics 



1. INTRODUCTION 

In many cosmic systems containing jets, e.g., ac- 
tive galactic nuclei (AGNs), stellar mass black hole X- 
ray binaries (BHXRBs), and young stellar objects, the 
presence of an accretion disk is inferred from observa- 
tions. One of the prominent theories of jet produc- 
tion asserts that the jet plasma is propelled by mag- 
netic field lines that thread the accretion disk (and 
perhaps, in AGNs and BHXRBs, the black hole) and 
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are t w isted by differential rotation iB landford fc Pavne 
1982 ; iLovelace. Berk fc Contopoulod .1991; Bcgelman 



1995; IMeier. Koide fc IJchidal 120011 : IVlahakis & Konig] 



2004f ). This suggests that any relations observed between 



events in the accretion disk and those in the jet can be 
used to characterize the disk-jet connection and there- 
fore to constrain models of jet formation. In BHXRBs, 
the connection between accretion state and events in 
the jet is well-established. In these objects, transitions 
to high-soft X-ray states are associated with the emer- 
gence of very bright fea tures that proceed to propa- 
gate down the radio iet ([Fender. Belloni fc Gallol 120041: 
iFender. Homan fc Bellonill2009[ ). 

Observations related to the disk-jet connection in 
AGNs are complicated, since a single AGN usually does 
not show the entire range of properties that we wish 
to analyze. For example, in Seyfert galaxies, most of 
the observed X-ray (and optical) emission is generated 
in the accretion disk-corona region, but their radio jets 
tend to be weak and non-relativistic fe.g.. lUlvestad et aLl 
119991 ). On the other hand, in radio- loud AGNs with 
strong, highly variable nonthermal radiation (blazars), 
the Doppler beamed emission from the jet at most wave- 
lengths masks the thermal emission from the accretion 
disk and its nearby regions. The radio galaxies 3C 120 
and 3C 111 provide an excellent opportunity to study 
the relationship between events in the accretion disk and 
those in the jet. At optical and X-ray frequencies, their 
properties are similar to Seyfert galaxies, e.g., the pres- 
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ence of strong broad eraission lines in the optical band 
(|Sargentlll977f). and an iron emission line in the X-ray 
spectrum ( Lewis et alj I2005D , while at centimeter and 
millimeter wavelengths their emission is dominated by 
a bright relativistic jet. Nevertheless, the source of the 
optical continuum emission in radio galaxies is not well- 
known. For exampl e, it may be bl ackbody emission from 
the accretion disk ()Malkanlll983() or synchrotron radia- 
tion from the jet. 

An accreti on disk-jet conn e ction in AGNs has been 
suggested by iMarscher et al.l (|2002l ). During '^3 yr of 
monitoring of the radio galaxy 3C 120 (1997 to 2000), 
each of four dips in the X-ray flux, accompanied by 
spectral hardening, preceded the appearance of a bright 
knot moving down the radio jet at a superluminal ap- 
parent speed. Using similar, but more extensive, multi- 
frequency monitoring observations of this object between 
2002 and 2007 at X-ray (2.4-10 keV) and radio (37 
GHz) wave bands as w ell as VLBA images at 43 GHz, 
IChatterjee et al.l ()2009( ) established an unambiguous con- 
nection between events in the accretion disk and jet. 
They showed that, during this interval, the X-ray and 
37 GHz variations are anti-correlated, with the X-ray 
variations leading those in the radio by ~120 days. Fur- 
thermore, nearly every X-ray dip is followed by the ejec- 
tion of a new knot in the VLBA images. These findings 
imply that a decrease in the X-ray production is linked 
with increased speed in the jet flow, causing a shock wave 
to form and propagate downstream, appearing as a knot 
of bright emission in the jet. In 3C 120, the presence 
of an X-ray Fe emission line implies that most of the 
X-rays are produced in the accretion disk region. There- 
fore, the above pattern is strong evidence for a disk-jet 
connection. 

In this work, we present the results of extensive multi- 
frequency monitoring of 3C 111 between 2004 and 2009 
at X-ray energies (2.4-10 keV), optical R band, and ra- 
dio frequencies 14.5, 37, and 230 GHz, as well as imag- 
ing with the Very Long Baseline Array (VLBA) at 43 
GHz. We use these data to investigate the presence of a 
disk-jet connection similar to that in 3C 120. 3C 111 is 
a relative ly nearby (z = 0.049) broad-line radio galaxy 
(BLRG) (jSargenti 11971 ILewis et al.l[2005h . It is clas- 
sified as a Fanaroff-Riley class II (FR II) radio galaxy 
|Fanaroff & Riley 1974) containing two radio lobe s wit h 
hot spots and a sing le-sidedjet (ilinfield fc Perleylll98l . 
At radio frequencies, 3C 111 has blazar-like behav- 
ior, although the jet lies at an angle ~18° to our 
fine of sight (jJorstad et al.ll2005l ). significantly greater 
than is the case for typical blazars. The counter-jet 
is too faint to detect, presumably because of Doppler 
de-boosting. 3C 111 exhibits the brightest compact 
radio core at cm/mm wavelengths of all FR II radio 
galaxies and a bla z ar-like spectral ene rgy distribution 
(|Kadler et all 120081: ISguera et al.ll2005[) . It was one of 
the first radio galaxies in which superluminal motion was 
detec ted (jGotz et al. 1987; Preuss, Alcf, & Keller mann 
119871 ): knots are ejected 1-2 times per year with typi- 
cal ap parent speeds of 3-5c (f Jorstad et al.l |2005() . Re- 
cently, iHartman. Kadler fc Tuelleri (|2008[ rshowed that 
part of the 7-ray emission from the EGRET source 
3EG J0416-H3650 is associated with 3C 111, and inte- 
gration over the first 11 months of observations with 
the Large Area Telescope of the Fermi Space Tele- 



scope detected 3C 111 at a rather low flux level, 1.5 ± 
0.5 X 10~^ phot cm ~^ s^^ at energies exceeding 1 GeV 
(jAbdo et all 120100 . At optical and X -ray frequencies, 
3C 111 possesses properties similar to Seyfert galaxies 
and BHXRBs. It has a promine nt iron emission line 
at a rest energy of 6.4 keV (e.g., iRevnolds et al.l[T998l : 
lEracleous. Sambruna fc Mushotzskvll2000( ) . This implies 
that most of the X-rays are produced in the immediate 
environs of the accr etion disk: the corona, a hot wind, or 
the base of the jet (iH aardt, Maraschi fc GhiseUini' 119941: 
iMarkoff. Nowak. fc Wilms, 2005: Navakshin.2000) . 

Power spectral density (PSD) analysis is a common 
technique to characterize time variability, which in turn 
is a powerful diagnostic of the geometry and physics of 
the accretion disk, corona, and jet. The PSD corresponds 
to the power in the variabili ty of emiss ion as a function 
of tim escale. iMarshall et al.l (2009) and lChatteriee et al.l 
([2009f l found that the X-ray PSD of 3C 120 can be fit 
by a piece-wise power law, with a slope that steepens 
above a break frequency. This property of 3C 120 is 
similar to Seyfert galaxies and BHXRBs, although in 
many cases the PSDs of the latter have more than one 
break, as well as sharp peaks representing quasi-periodic 
oscillations. The PSD break frequency in BHXRBs 
and Seyfert galaxies scales inversely with the m ass of 
the black hole (Edelson & Nandra 1999; Nowak ~t al.l 
1999'; Pounds et al. 2001; Uttley, McHardy fc Papadakij 
2002; Markowitz et al . 2003; McHardy et al. 2004^. 
McH ardv et al.l ()2006t ) showed that the break frequency 
and BH mass are inversely related over a range of BH 
masses from 10 to 10^ Mq, and that, for a given black 
hole mass, the break frequency also increases with higher 
accretion rate (m). 

Here, we use our data set to verify that the X-ray 
PSD of 3C 111 also has a break. Measuring the PSD 
of 3C111 gives us an opportunity to test whether an- 
other radio-loud AGN (in addition to 3C 120) has a 
broadband X-ray PSD similar to those of other Seyfert 
galaxies and BHXRBs. We can investigate whether the 
break timescale (Tb) is consistent with the empirical re- 
lation between Tb , black hole mass (AIb h). andbolomet - 
ric luminosity (Lboi) proposed bv lMcHardv et al.l (|2006f l. 
based on a sample of mainly radio-quiet Seyfert galaxies. 
We analyze the X-ray and radio light curves as well as 
the sequence of VLBA images to test whether the X-ray 
dips and ejection of new radio knots in the pc scale jet 
are temporally related, as is the case in 3C 120. In order 
to investigate the source of optical emission in 3C 111, 
we cross-correlate the optical flux variations with those 
at X-ray energies. Uneven sampling, as invariably oc- 
curs, can cause the correlation coefficients to be artifi- 
cially low. In addition, the time lags can vary across 
the years owing to physical changes in the source. In 
light of these issues, we estimate the significance of the 
derived correlation coefficients by repeating the analy- 
sis with simulated light curves, based on the underlying 
PSD. We also follow the variation in the degree of po- 
larization of the optical emission in order to discriminate 
between nonthermal and thermal emission. 

In §2 we present the observations and data reduction 
procedures, while in §3 we describe the correlated vari- 
ations in the X-ray continuum and Fe Ka emission line. 
In §4 we carry out the power spectral analysis and dis- 
cuss its results and implications. In §5 we investigate 
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Table 1 

Parameters of the Light Curves. 





Data set 


Start 


End 


T(days) ^ 


AT (days) ^ 


"NT C 




Monitor 


2004 March 


2010 April 


2246.0 


15.0 


822 


X-ray 


Medium 


2006 November 


2007 January 


56.0 


0.25 


224 




Longlook 


2009 February 16 


2009 February 17 


1.2 


0.01 


120 


Optical 


Monitor 


2004 November 


2010 August 


2080.0 


15.0 


336 




Medium 


2005 July 


2006 January 


160.0 


5.0 


37 


14.5 GHz 


Monitor 


2004 l-'ebruary 


2010 February 


2216.0 




398 


37 GHz 


Monitor 


2005 January 


2008 December 


1458.0 




98 


230 GHz 


Monitor 


2004 January 


2010 April 


2288.0 




161 



^ Total length of light curves. 
^ Bin size. 

Number of data points in the unbinned light curve. 
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Figure 1. X-ray light curves of 3C 111 with different sampling 
rates. The dotted and dashed vertical lines at the top panel denote 
the intervals shown in the middle and bottom panel, respectively. 



the relation between the X-ray dips and ejections of su- 
perluminal radio knots, as well as the X-ray/radio flux 
correlation. In §6 we cross-correlate the X-ray and opti- 
cal variations and discuss the possible sources of optical 
emission. In §7 we compare the properties of 3C 111 with 
those of the radio galaxy 3C 120, while §8 contains the 
summary and conclusions. 

2. OBSERVATIONS AND DATA ANALYSIS 

Table [T] summarizes the intervals of monitoring at dif- 
ferent frequencies for each of the three wave bands in 
our data set. We term the entire light curve "monitor 
data" ; shorter segments of more intense monitoring are 



described below. 

2.1. X-Ray Monitoring 

The X-ray light curves are based on observations of 
3C 111 with the PC A detector of the Rossi X-ray Tim- 
ing Explorer (RXTE) from 2004 March to 2010 April, 
with typical exposure times of 1-2 ks. For each ex- 
posure, we used routines from the X-ray data analy- 
sis software FTOOLS and the program XSPEC ver- 
sion 11.6 to calculate and subtract the standard X- 
ray background model for the faint sources from the 
data and to fit the X-ray spectrum from 2.4 to 10 keV 
as a power law with low-energy photoelectric absorp- 
tion by intervening gas in our Galaxy. For the lat- 
ter, we used a hydrogen column density of 9.0 x 10^^ 
atoms cm~^. Thi s value of t he absorbing column 
was d etermined by iSambruna. E racleous & Mushotzkvl 
()1999() by analyzing the X-ray spectrum of 3C 111 
with data from the Advanced Satellite for Cosmol- 
ogy and Astrop hysics fASC A ). Usi ng observations from 
XMM-Newton, ILewis et all (|2005[ ) have also obtained 
a value close to this. It should be noted that a gi- 
ant molec ular cloud is present in the foreground of 3C 
111 (.Marscher. Moore fc Banial[l99 3: Moo re k Marscheil 
Il995f ) and it contributes significantlv to the optical to X- 
ray extinction of emission received from 3C 111. 

The sampling of the X-ray flux varied. Normally, ob- 
servations were made 2-3 times per week except during 
8-week intervals each year when the radio galaxy is too 
close to the Sun's celestial position to observe safely. In 
order to sample shorter-term variations, between 2006 
November and 2007 January we obtained, on average, 
four measurements per day for almost two months. We 
refer to these observations as the "medium" data. Our 
monitoring program also contains an interval from 2008 
September to 2008 December when we obtained roughly 
one data point per day. In order to sample variations 
on as short timescales as possible, we observed 3C 111 
with the European Photon Imaging Camera (EPIC) on 
board XMM-Newton continuously for about 130 ks on 
2009 February 16-17. The EPIC-pn and EPIC-MOS 
data were obtained in the Small Window and Partial 
Window mode, respectively, using the thin filter. The 
data were processed with the the Science Analysis Soft- 
ware version 8.0.0. We filtered the pn data to include 
only single and double-pixel events (i.e., PATTERN < 
4). A pn light curve was extracted in the 2.4-10 keV 
energy band, similar to that of our RXTE data, and was 
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Figure 2. Light curves for 2.4-10 keV continuum flux F2- 



(in 



erg cm ^ s ^ ) , photon index Fx (of the continuum) , intensity of 
the Fe line /fcKq (in ph cm~^ and the Fe line equivalent 

width EQW (in eV) derived from time-resolved spectral fits to 
the RXTE-PCA monitoring data. The dashed lines are the mean 
values. 



background-subtracted and binned to 1000 s time inter- 
vals. In general, pn data are more suitable for timing 
analysis and the average count rate in the pn light curve 
was higher than that in the MOS light curves. Hence, we 
use the pn light curve as the "longlook" data. Figure [T] 
presents these three data sets. The X-ray data are given 
in Tables [2] and [3] of the electronic version. 

We performed time-resolved spectroscopy on the 
RXTE monitoring data in order to explore the vari- 
ability of the Fe Ka emission line. If the line does 
indeed respond to X-ray continuum variations, then 
one can place constraints on the distance between the 
X-ray continuum source and the bulk of the line- 
emitting gas. The analysis closely followed that em- 
ployed by previous studies (e.g., IV auglian &: EdelsonI 
120011: iMarkowitz. Edelson fc VaughanI 120031) . Summed 
spectra were constructed by combining the data from 
the individual '^1 ks observations within the defined 
time bins. A time-averaged model was applied to all 
binned spectra. The first five time bins were defined to 
match the start/stop times of the RXTE campaigns in 
1997 March, 1999 January-February, 2001 March, 2001 
September and 2003 September. For the sustained long- 
term monitoring that started on 2004 March 1, we con- 
sidered all data acquired up to 2010 February 14. Bin 
sizes were chosen to optimize the trade-off between min- 
imizing the uncertainties in the Fe line flux and maxi- 
mizing the number of bins. This yielded 31 bins each 
of roughly 60 days in duration. The response of the 



PCA and the conversion factor between count rate per 
PCU and incident flux, evolve slowly. In the 2-10 keV 
band, 1.0 ct/s/PCU corresponds to roughly 1.12x10"^^, 
1.09x10"", and 1.04x10"" erg cm-^ during 1997, 
2003, and 2008, respectively. Consequently, regarding 
long-term PCA light curves, fluxes obtained from spec- 
tral fitting are more valuable than count rates. We have 
taken this long-term trend in response into account in 
our analysis. The number of PCUs used during the anal- 
ysis of data from various epochs are as follows: PCU 0, 

1 and 2 for 1997 March, PCU and 2 for 1999 Jan-Feb, 
since PCU 1 had started to have thermal breakdown in 
1998, and PCU 2 only for all data from 2001 to 2010, 
since PCU lost its propane layer in 2000. The number 
of counts in a typical binned spectrum ranged from ^^10^ 
to 10^. 

We summed the PCA spectra from the individual ob- 
servations within each time interval to create the binned 
spectra. PCA response matrices were generated for 
each individual observation using pcarmf v. 11.7 and 
summed, weighted by exposure time, to yield a response 
for each binned spectrum. We did not include HEXTE 
data in the time-resolved fits, as those data yielded no 
additional model constraints compared to the PCA alone 
on these time scales given the flux level of 3C 111. Spec- 
tra were fit over the 3-25 keV energy range in XSPEC 
V. 12.5. Ik using a model of the same form as the best-fit 
to the long-term time-averaged spectru m determined by 
IRivers. Markowitz fc Rothschildl (j201lD . Those authors 
performed joint RXTE PCA + HEXTE spectral fitting 
over the 3-100 kcV energy band to a summed spectrum 
comprising all RXTE data on 3C 111 observed from 1997 
March until 2008 July. The model we use thus consists of 
a power law with a photon index Fx allowed to vary after 
being initially set at 1.67, plus a Gaussian component to 
model the Fe Ka emission line, with photo-electric ab- 
sorption in our Galaxy similar to that given above. The 
intensity of the Gaussian component /poKa was allowed 
to vary; the energy centroid and width a were kept fixed 
at 6.19 keV (rest-frame) and 0.72 keV, respectively, since 
there was no significant evidence to allow these param- 
eters to vary from these best-fi t time -averaged values. 
IRivers. Markowitz fc Rothschildl (1201 ID found no signif- 
icant evidence to warrant inclusion of a Compton hump 
in the model. Background level corrections were applied 
using RECORN within XSPEC. 

Figure [5] shows the resulting light curves for the 2-10 
keV continuum flux ^2-10, Fx, /fcKq, and the Fe line 
equivalent width EQW . These data are given in Table |4] 
of the electronic version. Errors on /fcKq, Fx, and EQW 
were derived using the point-to-point varianc 

We de- 
termined errors on values of i^2-io within a time bin from 
the 2-10 keV continuum light curve, using the mean flux 
error on the '^1 ks exposures in that time bin. 

2.2. Optical Monitoring 

We monitored 3C 111 in the optical R band over a 
portion of the time span of the X-ray observations. The 
majority of the measurements in R band are from the 

2 m Liverpool Telescope (LT) at La Palma, Canary Is- 



" See §3.3 of IVaughan fc EdelsonI 1 120011) or §3.3 of 
IMarkowitz et aLl II2003I ) for further details and the definition of 
the point-to-point variance. 
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Table 2 

X-ray (2.4-10 KeV; RXTE-PCA) light curve of 3C 111 from 2004 to 2010 (first 5 rows). 





exp (s) 


fiux'^ 


err 


a'' 


err 


counts/s/PCU 


err 


3065.9070 


2080. 


2.565E-11 


5.683E-13 


-0.691 


0.062 


2.613EH 


-00 


5.307E-02 


3068.8816 


752. 


1.975E-11 


8.709E-13 


-0.530 


0.123 


2.035EH 


-00 


8.308E-02 


3072.8645 


1952. 


1.503E-11 


5.269E-13 


-0.672 


0.096 


1.544EH 


-00 


4.893E-02 


3075.8169 


2032. 


1.458E-11 


5.173E-13 


-0.657 


0.100 


1.480EH 


-00 


4.770E-02 


3080.0244 


1936. 


1.292E-11 


5.291E-13 


-0.394 


0.112 


1.329EH 


-00 


4.920E-02 



Relative Julian Date: JD-2450000. 
Unit: ergs cm~^s~^. 
a: Energy Index. 



Table 3 

X-Ray Longlook light curve 
of 3C 111 from 
XMM-Newton in 2.4-10 
keV band on 2009 Feb 
16-17 (first 5 rows). 



RJD^' 


Counts/s 


error 


4878.50 


4.26 


0.12 


4878.51 


4.06 


0.11 


4878.52 


3.91 


0.12 


4878.53 


3.98 


0.13 


4878.55 


3.98 


0.13 



Relative Julian Date: 
JD-2450000. 



lands, Spain, supplemented by observations at the 1.8 m 
Perkins Telescope of Lowell Observatory, Flagstaff, Ari- 
zona. We used star A (Figure [S]) located ~50" west of 
3C 111 as a comparison star, for which we determined an 
i?-band magnitude R — 14.93 ±0.05 based on differential 
photometry with standard stars from other fields of view. 
Data acquisition and analysis procedures wer e identical 
to those described bv iChatteriee et all (|2009[ ). We also 
monitored the i?-band linear polarization of 3C 111 from 
2006 to 2009 with the PRISM camera on the Perkins 
telescope at Lowell Observatory and with the 2.2 m tele- 
scope at Calar Alto Observatory, SpairFl. Interstellar 
polarization (ISP) is significant in the field of 3C 111. 
We used the polarization parameters of the comparison 
star, presumed to be intrinsically unpolari zed, to correct 
the p olarization of 3C 111 for ISP (see also lJorstad et all 
|2007| ). We corrected the degree of li near polarization for 
statistical bias using the method of iWardle fc Kronberd 
^974) . 

2.3. Radio Monitoring 

We monitored the 14.5 GHz fiux of 3C 111 with the 26 
m antenna of the University of Michigan Radio Astron- 
omy Observatory. Details o f the calibration a nd analysis 
techniques are described in lAller et all (119850 . The flux 
scale is set by observations of Cassiopeia A ()Baars et all 
Il977f l. We have compiled a 37 GHz light curve with 
data from the 13.7 m telescope at Metsahovi Radio Ob- 
servatory, Finland. The flux density scale is based on 
observations of DR 21, with 3C 84 and 3C 274 used 
as secondary calibrators. A detailed description of the 

Data acquired as part of the MAPCAT program: 
[http: / /www.iaa.es/f^iagudo/resea rch/MAPCAT 




Figure 3. Finding chart for 3C 111. Star A is the comparison 
star with R = 14.93 ± 0.05, located --50" west of 3C 111, used 
to calculate the magnitude and interstellar polarization correction 
(p = 3.6 ± 0.2%, X = 123° ± 2°). 



data r eduction and analysis is given in iTerasranta et afl 
(|1998D . We have also compiled a 230 GHz light curve of 
3C 111 over a portion of the time span of the X-ray obser- 
vations with cali bration data from the the Submillimeter 
Array (SMA; see lGurweh et al.ll2007l for details). 

Figure [4] presents the X-ray, optical, and radio light 
curves. The variation of the total R-band flux, po- 
larization percentage, and electric vector position angle 
(EVPA) are shown in Fig. \5\ These data are given in 
Tables [Si El [3 HI and [9] of the electronic version. 

2.4. VLBA Imaging 

Starting in 2004 May, we observed 3C 111 with the 
Very Long Baseline Array (VLBA) at 43 GHz at roughly 
monthly intervals, with some gaps of 2-4 months. The 
sequence of images (Figure (5]) from these data provide 
a dynamic view of the jet at an angular resolution ~0.1 
milliarcseconds (mas) in the direction of the jet, corre- 
sponding to 0.094 pc for an adopted Hubble constant of 
Hq — 70 km s~^ Mpc~^. We processed the data in the 
same manner as described in iJorstad et al.l (|2005l ). We 
modeled the brightness distribution at each epoch with 
multiple components with circular Gaussian brightness 
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Table 4 

Variation of 2.4-10 keV continuum fiux, intensity of the Fe line, 
photon index, and the Fe Une equivalent width derived from 
time-resolved spectral fits to the RXTE-PCA monitoring data shown 
in Figure [2] 



Yr 


f a 
J^cont 


err 


f b 
tlinc 


err 




err 


EQW^ 


err 


1997.22 


4.26 


0.03 


6.05 


2.39 


1.71 


0.06 


104.68 


41.40 


1999.05 


2.25 


0.09 


2.12 


2.39 


1.59 


0.06 


67.64 


76.44 


2001.20 


7.47 


0.07 


9.14 


2.39 


1.77 


0.06 


90.74 


23.76 


2002.69 


5.00 


0.04 


7.24 


2.39 


1.72 


0.06 


106.76 


35.29 


2003.67 


5.11 


0.03 


10.13 


2.39 


1.77 


0.06 


148.30 


35.03 



2-10 keV continuum flux F2-10 (10 erg cm ^ s ^). 
^ Intensity of the Fe line /peKa (10~^ph cm"^ 
photon index Fx of the continuum. 
Fe line equivalent width EQW (eV). 



Table 5 

Optical (R Band) light curve of 3C 111 from 2004 to 
2010 (first 5 rows). 



Table 8 

Radio (230 GHz) light curve of 3C 
111 from 2004 to 2010 (first 5 rows). 



RJD^' 


R Mag 


err 


Flux Density'' 


err 


3325.5154 


17.033 


0.009 


0.474 


0.004 


3325.5171 


17.090 


0.009 


0.450 


0.004 


3325.5210 


17.047 


0.009 


0.468 


0.004 


3325.5247 


17.055 


0.009 


0.465 


0.004 


3325.5266 


17.030 


0.009 


0.475 


0.004 



Relative Julian Date: JD-2450000. 
^ Unit: mjy. 



RJD'' 


Flux Density (Jy) 


error 


3022.48 


1.26 


0.09 


3031.60 


1.71 


0.35 


3266.67 


2.06 


0.13 


3359.38 


2.49 


0.13 


3381.64 


3.80 


0.20 



Relative Julian Date: JD-2450000. 



Table 6 

Radio (14.5 GHz) light curve of 3C 
111 from 2004 to 2010 (first 5 rows). 



Table 9 

Variation of the R-band Polarization 
Percentage and EVPA of 3C 111 from 2006 
to 2010. 



RJD'' 


Flux Density (Jy) 


error 


RJD^' 


P (%)'^ 


error 


EVPA'^ 


error 


3024.5 


2.97 


0.10 


3741.74 


3.32 


0.60 


42.70 


5.20 


3042.5 


3.06 


0.04 


3796.65 


0.60 


0.80 


97.50 


38.30 


3050.5 


3.16 


0.04 


4048.87 


3.33 


0.70 


65.80 


6.00 


3052.5 


3.09 


0.04 


4049.98 


3.06 


0.40 


22.70 


3.70 


3064.5 


3.02 


0.03 


4076.88 


1.49 


0.80 


26.00 


15.40 



Relative Julian Date: JD-2450000. 



Relative JuUan Date: JD-2450000. 
^ R-band polarization percentage. 



Electric vector position angle. 



Table 7 

Radio (37 GHz) light curve of 3C 
111 from 2005 to 2008 (first 5 rows). 



RJD^ 


Flux Density (Jy) 


error 


3371.22 


4.61 


0.25 


3372.20 


4.35 


0.24 


3379.52 


4.19 


0.22 


3423.03 


4.23 


0.23 


3447.93 


4.17 


0.26 



Relative Julian Date: JD-2450000. 

distributions usi ng the task MODELFI T of the software 
package Difmap (jShepherd et al.l[l997l ). This represents 
the jet emission downstream of the core, the bright fea- 
ture situated at the western end of the jet, by a sequence 



of knots (also referred to as "components"), each charac- 
terized by its flux density, FWHM diameter, and position 
relative to the core. The ejection time Tq is the extrapo- 
lated time of coincidence of a moving knot with the posi- 
tion of the (presumed stationary) core. Figure[7]plots the 
distance vs. epoch for all moving components brighter 
than 100 mJy within 2.0 mas of the core. We use the po- 
sition vs. time data to determine the projected direction 
on the sky of the inner jet, as well as the apparent speeds 
and ejection times of new superluminal knots. Table [TU] 
lists the values of Tq and apparent speeds of new knots 
determined by the above procedure. In the top panel of 
Figure [H the arrows represent the times of superluminal 
ejections, while the line segments perpendicular to the 
arrows show the uncertainties in the values of Tq. The 
apparent speeds of the components with well-determined 
motions are between 2.2 and 6.4c. 
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Figure 4. Variation of X-ray flux, optical flux density, and radio flux density of 3C 111 from 2004 to 2010. In the top panel, the 
arrows indicate the times of superluminal ejections and the line segments perpendicular to the arrows represent the uncertainties in the times. 



Table 10 

Times of X-ray Dips, Ejection Times (To) of Radio Knots, 
and Apparent Speeds (/3app) of Radio Knots. 



Dip ID 


Txmin" 


Knot ID 


To 


/?app 


dl 


2004.19 


Kl 


2004.30± 0.14 


3.99 ± 0.19 


d2 


2005.02 


K2 


2005. 18± 0.05 


6.37 ± 0.29 


d3 


2005.57 


K3 


2005. 70± 0.04 


3.44 ± 0.08 


d4 


2006.47 


K4 


2006.61± 0.09 


2.24 ± 0.44 


d5 


2006.88 


K5 


2007.01± 0.04 


2.71 ± 0.38 


d6 


2007.32 


K6 


2007.66± 0.09 


4.18 ± 0.11 


d7 


2008.51 


K7 


2008.83± 0.07 


4.54 ± 0.38 


d8 


2008.98 


K8 


2009.07± 0.08 


4.07 ± 0.43 


d9 


2009.26 


K9 


2009.29± 0.04 


4.33 ± 0.66 



^ Time of X-ray minimum. 



3. CORRELATED VARIATIONS IN THE X-RAY 
CONTINUUM AND EE Ka EMISSION LINE 

Figure |8] shows the zero-lag correlation diagrams of 
-^FcKq, Tx, and EQW, each plotted against F2_io- Best- 
fit linear relations are plotted for each and are summa- 
rized in Table [TT] along with the Pearson correlation coef- 
ficients. From this diagram, we can see that the range of 
J^2-io is from 2x10"" to 8x10"" (factor of -4), while 
that of IpcKa is from 4xlO~^ to 12x10"^ (factor of ~3). 



Hence, we conclude that at least 75% of the Fe line flux 
responds to X-ray continuum variations. We searched 
for lags between F2-10 and /pcKa w ith the interpolated 
cross correlati on function (IC CF; iGaskell fc Peters"onl 
Il987t [white fc PetersonI 119941) and the discrete cross - 
correlation function (DCCF; lEdelson fc KroUkI [1988), 
with errors determined using the Flux-Randomization 
and Ran dom Subset Selectio n (FR-RSS) technique pro- 
posed by [Peterson et all (|1998 ). We omitted the first five 
points (1997-2003) from this analysis due to the large 
data gaps present. The ICCF peak correlation coeffi- 
cient was 0.61 at a continuum-to-line delay r = days, 
with a 1(7 upper limit to the lag of |r| — 88 days (Figure 
[9]). This implies that the bulk of the line flux originates 
within ~90 light-days of the X-ray continuum source if 
it lies near the plane of the accretion disk. Figure [8| also 
demonstrates that Fx increases with i^2-io, similar to the 
behavior seen in many Seyfert galaxies t hat lack a strong 
jet contribution to the X-ray spectru m ([Papadakis et al.l 
[200l:IShih. Iwasawa fc Fabianll2002D . 

The X-ray continuum/Fe line correlation presented 
here is perhaps the strongest such correlation observed 
so far (cf., for example, t he continuum-line corr elation 
observed in NGC 3227 bv [Markowitz et al.[[200l where 
'--^50% of the line flux responded to the continuum 
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Figure 5. Variation of the optical R-band flux density, polar- 
ization percentage, and electric vector position angle (EVPA) of 
3C 111 during 2006 to 2010. The average polarization was only 
1.6% ± 0.6%, rather low for synchrotron radiation from the jet. 
This supports the notion that most of the optical emission is 
thermal radiation from the accretion disk. 



Table 11 

Best-fit Linear Relations For Parameters Derived From 
Time-resolved Spectroscopy. 



4. POWER SPECTRAL ANALYSIS 
4.L X-ray 

We use a varia nt of the Power Spectrum Response 
method fPSRESP: lUttlev. McHardv k Papadakisll200l) 
to determine the intrinsic PSD of the X-ray hght 
curve. Our re a lizatio n of PSRESP is described in 
iChatteriee et all (|2008f ). PSRESP gives both the best- 
fit PSD model and a "success fraction" Fgucc (fraction 
of simulated light curves that successfully represent the 
observed light curve) that indicates the goodness of fit of 
the model. 

At first we fit a simple power-law model to the X-ray 
PSD, but found that the value of i^succ is unacceptably 
low (0.1). This implies that a simple power law is not 
the best model for this PSD. We therefore fit a bend- 
ing power-law model (broken po wer law with a smooth 
break) to the X-ray PSD (see iMcHardv et all 120041: 
IChatteriee et al.ll2009[) . 



Note. — The parameters for the best-fit linear solutions plotted in 
Figure [S] rp is the Pearson correlation coefficient and Pr is the null 
hypothesis probability, i.e., the probability of obtaining the observed 
correlation by chance. 



variations on time scales of < 700 days). If the bulk 
of the Fe line originates within 90 light-days of the 
X-ray continuum source, then an origin in material com- 
mensurate with the optical broad line region (located 
light-days to light-weeks away from the central black 
hole in Seyferts) is plausible. However, this delay is 
inconsistent with a model wherein the bulk of the line 
emission originates in a '^parsec-scale, homogeneous 
molecular torus commonly invoked in Seyfert 1/2 
unification schemes. 



(1) 



Parameters 


Slope 


y-intcrcept 




Pr 


iFeTia VS. ^2-10 


1.27 ± 0.31 


1.83 ± 1.65 


0.594 


1.3 X 10-* 


Px vs. F2_10 


0.040 ± 0.008 


1.52 ± 0.04 


0.796 


5.8 X 10"* 


EQW vs. F2-10 


-5.51 ± 4.97 


150 ± 29 


-1.92 


0.26 



Here, A is a normalization constant, i>b is the break 
frequency, and aH and are the slopes of the power 
laws above and below the break frequency, respectively. 
During the fitting, we varied the break frequency vb 
from 10"^ to 10"^ Hz in steps of 10°"^, an (slope of 
the power law above the break frequency) from —1.5 
to —3.0 in steps of 0.1, and ckL (slope below the break 
frequency) from —1.0 to —1.5 in steps of 0.1. These 
ranges include the values of a found in the light curves 
of BHXRBs, for which qt. ~ —1 and qh is between —2 
and -3 re.g.. lRemillard fc McClintoci3[2006h . This pro- 
cedure yields a much higher success fraction than the 
simple power-law model. Based on the model with the 
highest success fraction (0.85), we obtain a good fit with 
the parameters ql = — 1.0 ± 0.1, an — and 
logio(^'B) = — 6.05^Q3Q Hz, which is equivalent to a 
timescale 13^g^ days. Figure [TU] presents this best-fit 
model and the corresponding PSD. As seen in the fig- 
ure, the high frequency part of the PSD is dominated by 
Poisson noise. That is because i) this part of the PSD 
is generated from the longlook light curve, and fiuxes in 
the longlook light curve have larger uncertainties owing 
to shorter exposure times than those in the other light 
curves and ii) due to the red noise nature of the PSD, the 
intrinsic power is smaller at the higher frequencies and 
hence the power generated by Poisson noise is relatively 
more important. The figure shows that when the esti- 
mated Poisson noise is added to the best-fit model PSD, 
it matches the observed PSD quite well. 

The PSD break frequency in BHXRBs and Seyfert 
galaxies scales inversely with th e mass of the black hole 
(lUttley. McHardv fc PapadakisI 120021: IMcHardv et al " 
120041. I2n0fil: lEdelson & Nandral llQfla iMarkowitz et al 



2003). McHardy ct al. (2006) showed that this holds over 
a range of BH masses from 10 to 10* Mq, and that, for a 
given black hole mass, the break frequency increases with 
accretion rate (m). The bolometric luminosity (Lboi) can 
serve as a proxy for to. Since the BH mass in 3C 111 is 
not well-established, we need to derive its value from 
independent observations. To do this, we need to deter- 
mine the optical luminosity of the central engine. 
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Figure 6. Sequences of VLBA images at 7 mm during 2004 to 2009. Tlie images are convolved witti an elliptical Gaussian beam of 
FWHM size 0.32x0.16 mas at PA = -10°. The global peak over all maps is 4.80 Jy/Beam. The contour levels are 0.25, 0.354, 0.5, 0.707, 
90.51 % of the global peak. Note that the time spacing of the images is not uniform, hence the lines denoting the proper motion of a 
given knot do not pass through the centroid of the knot at every epoch. 



The optical emission received from 3C 111 is sub- 
ject to substantial extinction owing to the presence of 
a translucent (^v between ~1 and 5) molecular cloud 
in the foreground. The standard method for deter- 
mining the extinction, star counts in the surround- 
ing region on the sky, is subject to systematic uncer- 
tainties because of possible gradients in dust ab sorp- 
tion. In fact, iMarscher. Moore fc Banial ()1993( ) and 
IMoore fc Marscherl (jl995[ ) have reported time variations 
in radio H2CO absorption lines toward 3C 111, indicat- 
ing that the column density and / or excitation conditions 



may vary throughout the cloud, causing changes in ab- 
sorption as our line of sight to the quasar drifts across 
the cloud fro m rel a tive m otion between the Earth and 
the cloud. Sar gentI (jl977[ ) gives essentially simultaneous 
multi-band optical and near-infrared flux measurements 
of 3C 111, which are not corrected for extinction. We an- 
alyze this spectrum to calculate a value of the extinction, 
Av = 2.5, that generates an extinction-corrected spec- 
trum that can be fit well by a power law, the functional 
form that fits the optical continuum of a typical AGN. 
The spectral index of the extinction-corrected spectrum, 
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Figure 7. Angular separation from the core vs. epoch of the 
moving knots brighter than 100 mjy within 2.0 mas of the c ore. 
The black lines indicate the motion of each knot listed in Table [TOl 
A knot is identified through continuity of the trajectory from one 
epoch to the next. The identification of the knot across epochs 
considers its position angle relative to the core and, in some cases, 
its polarization. At some epochs, the model fit breaks a knot into 
two pieces. 



displayed in Figure [TTl is —1.38. Our value of is 
consistent with the values E(B-V) ~1 for translucent 
cloud s with total N(H) ^OxlO^i cm-^ (iRachford et al.l 
I2002D and Av/E(B- V) ^2.1 derived for the translucen t 
cloud HD 36982 by iLarson. Whittet fc Houghl (fToM ). 
We note, however, that these values have high uncer- 
tainties and could differ significantly from one cloud to 
another . Afte r correcting the optical luminosity from 
iSargenti (|1977[ ) by this value of A^, we arrive at a de- 
reddened optical luminosity of XLx = 6.1 x 10"'^ erg s~^ 
at A = 510 nm. 

In a recent paper. [Pecarli. Dotti fc Treves! (|2011[ ) used 
the properties of the Ha line of a sample of blazars in ad- 
dition to other parameters to determine the mass of the 
central BH. Using the deprojection factor and Ha line to 
continuum (510 nm) luminosity ratio from equation (3) 
and (5), respectively, of that paper , as well as the width 
and lu minosity of the Ha line from lEracleous &: Halpernl 
(j2003l) . we calculate the BH mass of 3C 111, Mbh = 
2.4tQ;5 X 10^ Mq. Another method for estimating 
the BH mass is to use the relationship between A/bh 
and the FWHM line width of the H/9 broad emission 
line along with equation (2) of IVestergaard k, PetersonI 
([2006). Unfortunately, no accurate measurement of 
FWHM(H^) in 3C 111 has been pubhshed; instead, we 
assume that FWHM (H^) » FWHM(Ha) 4800 km s'^ 
(|Eracl eous fc Halpern 2003). From this value and the 
de-reddened optical luminosity given above, we derive 
ATbh = 1-5^0^3 X 10* M0. If we use the more conservative 
estimate of the uncertainty in the zero point value of the 
scaling relation between the masses determined using re- 
verberation m apping and those from single-epoch spectra 
from iVestergaard fc Peterson. (,2006il . then the respective 
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Figure 8. Zero-lag correlation diagrams for intensity of the Fe line 
IycKc I photon index Fx (of the continuum) , and Fe line equivalent 
width EQW each plotted against 2-10 keV continuum flux F2—1Q 
with the best-fit linear relation plotted for each as a dashed line; 
units are the same as in Figure [2] The X-ray continuum/Fe line 
correlation presented here is one of the strongest such correlation 
observed so far. The best-fit linear relations plotted here are sum- 
marized in Tablc fm along with the Pearson correlation coefficients. 
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Figure 9. Points denote discrete cross-correlation function 
(DCCF), and the solid line shows the interpolated cross- 
correlation function (ICCF) of the variation of intensity of the 
Fe line /pcKa with that of 2-10 keV continuum flux F2^\o. The 
two variations are strongly correlated with a time delay consistent 
with zero with a la upper limit of 88 days. 
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Figure 10. Result of application of the PSRESP method to 
the X-ray light curve of 3C 111. The PSD of the observed 
data at high, medium, and low frequencies is given by the solid, 
dashed, and dotted jagged lines, respectively, while the underlying 
power-lawr model is given by the thicker solid bent line. Points 
with error bars (open squares, filled circles, and open circles for 
high, medium, and low frequency range, respectively) correspond 
to the mean value of the PSD simulated from the underlying 
power- law model (see text). The error bars are the standard 
deviations of the distribution of simulated PSDs. The broadband 
power spectral density is best described by a bending power law 
with low frequency slope —1.0, high frequency slope —2.5, and 
break frequency lO"^ "^ Hz. The high frequency part of the 
PSD is dominated by Poisson noise. The figure shows that when 
the estimated Poisson noise is added to the underlying bend- 
ing power-law model PSD, it matches the observed PSD quite well. 
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Figure 11. The extinction corrected optical spectrum of 3C 111. 
The best-fit power-law has a slope —1.38 for Av = 2.5. 



uncertainty in the value of the BH mass becomes larger, 
-^BH = l-5ti;o ^ ^o- Using the correlation between 
the wi dths of Ha and H/j emission lines from equation 
(3) of IGreene fc Hoi (pOOSi) . we calculate FWHM(H^) 
of 3C 111 to be 5400 ± 400 km/s from F WHM(Ha) = 
4800 km s-i (jEracleous fc Halpern| [ 2003D . Using this 
value in equation (2) of Vcstergaar d fc PetersonI (j2OO60 
and keeping other parameters same as above, we obtain 
Mbh ^ l.Stj^j X 10^ Mq. 

iMarchesini. Celotti fc Ferraresd (|2004[ ) derived a much 
higher value of Mbh estimating the bolometric luminos- 
ity of 3C 111 from the optical nucle ar luminosity alon g 
with the bolometric correction from lElvis et al.l ()1994[ ). 
The former is derived by fitting the unresolved nuclear 
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Figure 12. The correlation of H/3 line-width with the PSD 
break timescale for a sample of mainly radio-quiet AGNs from 
IMcHardv et al.| II2006I ) is denoted by the dashed straight line 
and the solid circles. The open circle and the open square show 
data points corresponding to radio galaxies 3C 120 and 3C 111, 
respectively, which follow the correlation. For 3C 111, line width 
of H/3 has be e n calc ulated from that of Ho using equation (3) of 
IGreene fc ilcl ||2005|) . 



component with the appropriate HST PSF in the re- 
spective HS T image. They adopted an e x tincti on cor- 
rection from ISchlegel. Finkbeiner fc DavisI (jl998[) which 
is higher than the extinction we calculated above by 
almost 3 ma gnitudes. If we adjust the calculation of 
iMarchesini. Celotti k Ferraresd ()2004D to our extinction 
correction of Av=2.5, we obtain Lboi = (2.5 ± 1.0) x 10^^ 
ergs s-^ We also calculat e Lboi = ( 4 8 ± 3 .3) x 10''^ 
ergs s^^ from Young, Elvis, fc Risalitil (|2010[ ) using the 
2.4-10 keV X-ray flux and 510 nm optical flux, which 
is consistent with the above value within uncertainties. 
Performing a simi lar calculation using equation (11) of 
iLusso et all (pOlOl ). we find Lboi = (1.8 ±1.1) x 10^^ ergs 
s"""^, which is close to the above values as well. 
Using the intermediate values from above, Lboi = 



j-t-0.5 
-0.4 



X 10« M 



0; 



(4.8 ± 3.3) X lO^-* erg/s, and Mbh = 1- 
we find Eddington ratio LboiZ-'^cdd = 0.02 ± 0.01. We 
adopt the best-fit values and uncertainti es in the Tb— 
Mbh — Lboi relation proposed by McHar dv et al.l (j2006D , 
and the above values of Mbh and Lboi- We then calcu- 
late the expected break frequency to be jj^. If 
we use the more conservative estimate of the uncertainty 
in BH mass from above, we obtain l0-6-6±i i Hz for the 
same. Both of these are consistent with the value we have 
obtained from the best-fit PSD, logio(t^B) = -6.05lo Jq 
Hz, within uncertainties. This demonstrates that 3C 111 
also follows the Tb — Mbh — Lboi relation, the first FR H 
radio gala xy for w hich this has been tested. 

McHardiLeTaD p006) also showed that FWHM(H/3) 
in AGNs is strongly correlated with the observed PSD 
break timescale (Lb) given by 

log(Ts) = 4.2 \og{[FWHM{Hl3)]) - 14.43 (2) 

U sing the results of t he ab ove analysis and that contained 
in iChatteriee et aP (|2009l ). we can now add data points 
correspondin g to the radio galaxies 3C 120 and 3C 111 
in Figure 3 of lMcHardv et al.l (|2006[ ). resulting in Figure 
[T2l The plot shows that these two radio galaxies also 
follow the empirical Tb versus FWHM(H/3) correlation. 
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Figure 13. Result of application of the PSRESP method to the 
optical light curve of 3C 111. The PSD of the observed data at 
high and medium frequencies is given by the dashed and solid 
jagged lines, respectively, while the underlying power-law model 
is given by the solid straight line. Points with error bars (filled 
circles and open squares for high and medium frequency range, 
respectively) correspond to the mean value of the PSD simulated 
from the underlying power-law model (see text). The error bars 
are the standard deviations of the distribution of simulated PSDs. 
The broadband power spectral density is best described by a 
simple power law with slope —2.0. The smoothly bending dotted 
line is the underlying be st-fit bending power-law model of the 
X-ray PSD from F'ig. 1101 This shows that the X-ray variability 
has more power than that in the optical on shorter timescales 
while on longer timescales variability at the two wave bands tend 
to be similar. 

4.2. Optical 

We use the same method as above to calculate the 
PSD of the optical variability. The optical i?-band PSD 
of 3C 111 f Figure [T^ shows red noise behavior, i.e., there 
is higher amplitude variability on longer than on shorter 
timescales. Based on the model with the highest success 
fraction, the optical PSD is best fit with a simple power 
law of slope —2.0^^'^, for which the success fraction is 
0.52. During the fitting, we varied the slope from —1.0 
to —3.0 in steps of 0.1. The rejection confidence, equal 
to one minus the success fraction, is much less than 0.9, 
hence the model provides an acceptable fit to the PSD. 
We also fit a broken power law to the optical PSD, set- 
ting the low- frequency slope at —1.0 and allowing the 
break frequency and the slope above the break to vary 
over a wide range of parameters (10~* to 10~^ Hz and 
— 1.0 to —2.5, respectively) while calculating the success 
fractions. This gives lower success fractions than the sim- 
ple power-law model across the entire parameter space. 
Since we do not have any significantly long segment of 
the i?-band data when the sampling frequency was more 
than 2-3 days per week, we determine the optical PSD up 
to the highest variational frequency that can be achieved 
with the existing data. A better constraint on the ex- 
istence of a break in the optical PSD could be achieved 
with a broader range of sampled frequencies. 

The slope of the optical PSD in the range 10~^'^ to 
-[^Q-5.8 is significantly steeper than that in X-rays over 
the same range (~1.0). Hence, the ratio of shorter to 
longer timescale variability amplitude is smaller at opti- 
cal than at X-ray wavelengths. 

4.3. Excess Variance 



We calculate the fractional root mean squared variabil- 
ity amplitude, i.e., the "excess variance" normalized by 
the square of the mean fiux, as a measure of the variabil- 
ity in the X-ray and optical data. It is defined as 



where S is the variance, CTorr is the obse rvational uncer- 
tainty, and X is the mean of the da ta (Nan dra et al.llT997l : 
lEdelson et all [2001 IVaughan et al. 2003) . We find that 
the values of Fvar of the X-ray and optical variability for 
the entire data set are 0.28 and 0.14, respectively. This 
supports the above conclusion from the PSD analysis 
that the optical is less strongly variable than the X-ray 
flux at the longest timescales we probe. It also implies 
that optical variations cannot be completely responsible 
for driving all of the X-ray variability observed, if we 
assume that most of the X-ray emission we observe is 
not relativistically beamed. 

5. X-RAY/RADIO CORRELATION 

We can see in Fig. |4] that, from the middle of 2009 to 
the present, the X-ray flux undergoes only minor fluctu- 
ations about a mean value of 5 x 10^^^ erg cm~^ s^^. 
During this interval, the 230 GHz flux is at a very low 
level, while the 15 GHz flux is decreasing. Therefore, the 
high-state X-ray plateau coincides with a low state of the 
radio jet. We can extend this X-ray/radio connection to 
the observations of earlier years as well. The times of 
"ejection" of new radio knots in the jet (shown by the 
arrows) are related to the low X-ray states. We convolve 
the X-ray light curve with a Gaussian smoothing func- 
tion with a 10-day FWHM smoothing time to identify the 
major long-term trends. Figure [T3] shows the smoothed 
X-ray light curve and denotes the times of superluminal 
ejections with arrows. Based on inspection of this light 
curve, we consider an X-ray fluctuation to be a signiflcant 
dip if the smoothed X-ray flux falls by more than 30% 
and remains at a low level for longer than one month. 
We calculate the central time of each dip by determining 
the local minimum of the X-ray flux. We note that the 
minimum amplitude measured from unsmoothed data of 
every dip is below the average of the X-ray flux over the 
length of the light curve (4 x 10~^^ erg cm~^ sec~^). It 
can be seen that each of the ejections is preceded by a 
signiflcant dip in the X-ray flux. The prolonged low level 
of X-ray flux in 2005 is associated with knot K3, which 
was extended along the jet axis, as is apparent in the 
images starting in late 2006 (Fig. This suggests that 
a prolonged disturbance near the base of the jet resulted 
in an elongated superluminal knot. Hence, d3 and the 
two signiflcant dips just before and after that are consid- 
ered to be one dip with multiple branches and is related 
to K3. Detection of the ejection of a new knot requires 
the analysis of a sequence of VLBA images of sufficient 
duration to define the trajectory and speed of a knot, as 
is evident from Figures |6] and [T] No new bright, mov- 
ing knots after K9 appeared in subsequent imaging with 
the VLBA during the first 8 months of 2010 (images not 
shown here). 

The time delay between the minimum of the X-ray dips 
and the time of ejection of the corresponding superlumi- 
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Figure 14. X-ray light curve of 3C 111 as in Figure |4l with data indicated by small circles. The curve corresponds to the same 
data, smoothed with a Gaussian function with a 10-day FWHM smoothing time. The arrows indicate the times of superluminal 
ejections and t he l ine segments perpendicular to the arrows represent the uncertainties in the times. X-ray dips and VLBA knots 
listed in Table 1101 are marked with the respective ID numbers. Over the six years of observation, significant dips in the X-ray 
light curve are followed by ejections of bright superluminal knots in the VLBA images. This shows a clear connection between 
the radiative state near the black hole, where the X-rays are produced, and events in the jet. Dip d3 and the two significant dips 
just before and after that are considered to be one dip with multiple branches and is related to K3 which was also extended along the jet axis. 



nal knot is distributed between 0.03 and 0.34 yr (Table 
[T0|) with a mean of 0.15 ± 0.08 yr. We plot the times 
of ejection of new knots along with the corresponding 
times of X-ray dips in Figure [151 A straight line fits 
the data extremely well with small scatter, which indi- 
cates that there is a clear association between X-ray dips 
and superluminal ejections. The best-fit line through 
the points has a slope of 1 and y-intercept of 0.15, which 
is identical to the mean delay of 0.15 ± 0.08 yr given 
above. We have performed a numerical simulation to 
calculate the significance of the dip-ejection correlation. 
If we keep the times of dips fixed and choose one million 
sets of nine ejections, times of which are drawn from a 
uniform random distribution between 2004.3 and 2010.3 
(the duration of our observations), the probability that 
there is at least one ejection event following all 9 dips 
by 0.03 to 0.34 yr is 0.007. In a similar simulation, we 
found that the probability that all 9 dips will have ex- 
actly one following ejection is 3 x 10~^. This shows that 
the dip-ejection association in 3C 111 is highly signifi- 
cant. To further investigate the relationship between the 
characteristics of the X-ray dips and VLBA knots, we 
calculate the average flux over all epochs of the latter. 
We do this via modeling of the brightness distribution at 
each epoch with multiple components characterized by 
circular Gaussian brightness distributions, as described 
in §2.4. This shows that the average flux of the VLBA 
knots corresponding to the shallower dips (d2, d4, and d8 
in Fig. [T4|l is significantly smaller (0.3±0.1 Jy) than that 
for the larger dips (dl, d3, d5, and d7), 0.9 ±0.4 Jy. Dip 
d6 is also shallow, but this may be due to the absence of 
data during the Sun-avoidance gap, while d9 is of inter- 
mediate depth. Hence K6 and K9 were not included in 
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Figure 15. Times of ejec tion of new VLBA knots versus times of 
X-ray dips from Table [TOl The dashed line is the best-fit straight 
line through the points. It fits the data extremely well with small 
scatter. This indicates that there is a clear association between 
X-ray dips and superluminal ejections. The y-intercept of this line 
minus 2004.0 indicates the value of mean time delay between the 
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the average calculation. That the brighter knots corre- 
spond to more pronounced dips supports the conclusion 
that the X-ray minima and the ejection of knots in the 
jet are physically related. 

As determined b y the discrete cros s-correlation func- 
tion (DCCF; Edel son fc KrolikI ITosl ) shown in the top 
panel of Figure 1161 the X-ray flux variations are weakly 
correlated with those at 15 GHz in 3C 111. The peak 
X-ray versus 15 GHz DCCF is 0.35. We simulate X- 
ray light curves from the underlying PSD as determined 



14 



1 


X-rays vs. 14.5 GHz 




0.5 











**** +*+*♦. 


0.5 






1 


t \ 


x-rays vs. 37 GHz ' 


0.5 

0.5 


- vV/v 

ft 

+ 





1 


'■J 1 


x-rays vs. 230 GHz' 


0.5 











tit' *t,( 


0.5 







-300 -200 -100 100 200 300 

Time Delay (Days) 



Figure 16. Discrete cross-correlation function (DCCF) of the 
X-ray and radio monitoring data of 3C 111. The time delay 
is defined as positive if the X-ray variations lag those at radio 
frequencies. Top panel shows the correlation function for X-ray 
versus 15 GHz variations while the middle and bottom panels 
show the same for X-ray versus 37 and 230 GHz, respectively. The 
jagged dotted and dashed lines at the top of each panel denote 99 
and 95 percent significance levels of correlation, respectively (see 
text). Those at the bottom of each panel show the same for anti- 
correlation. This shows that the X-ray and 37 GHz variations are 
correlated while the other two pairs are not significantly correlated. 

in §3 and resample the simulated light curves with the 
sampling window of the actual X-ray light curve. Then 
we cross-correlate the resampled simulated X-ray light 
curves with the real 15 GHz light curve. The solid and 
dotted jagged lines in Figure [TBI indicate 99% and 95% 
extremes of the distribution of correlation coefficients as 
a function of time delay. 

The DCCF indicates that the X-ray flux variations are 
correlated with those at 37 GHz in 3C 111, with a peak 
value of 0.63 (middle panel of Figure [TO]) . The position of 
the peak of the correlation function indicates the relative 
time delay between the variations at the two wavelengths. 
The time lag of the peak corresponds to X-ray variations 
leading those at 3 7 GHz by 40l^ g days. We use the 
FR-RSS technique (jPeterson et al.|[l99 8) to calculate the 
mean value and the uncertainty of the cross-correlation 
time lag. The bottom panel of Figure [16] shows that 
the X-ray and 230 GHz variations are weakly correlated, 
with a peak DCCF of 0.30. 

We note that if the radio flux increases within tens of 
days after an X-ray minimum and the X-ray flux recovers 
from the minimum to increase as well during this time, 



then a relatively weak X-ray/radio correlation will result 
as seen in the data presented in this paper. The signif- 
icant correlation between the dips and ejections of new 
knots in the jet implies that increased activity in the ra- 
dio jet of 3C 111 is related to a temporary decrease in 
X-ray production. Similarly, the association between the 
high-state X-ray plateau with a low state of the radio jet 
after the middle of 2009 suggests that a persistently high 
X-ray flux is associated with a low radio state. 

It is possible that the "corona," where the X-ray emis- 
sion seen in AGN is thought to arise from Compton up- 
scattering of so fter accretion-disk photons, migh t be the 
base of the jet (|Markoff. Nowak. fc Wilmsl 120051 ). If this 
is the case, then the X-ray flux will be related to the 
number of electrons residing there and available for scat- 
tering to create X-rays. If the same number of electrons 
is injected into the base of the jet per unit time, then 
faster flow velocities will correspond to lower densities 
in the scattering region. The mass loading of the jet 
should also afltect the asymptotic Lorentz factor of the 
flow downstream if the jet is magnetically driven (e.g., 
iVlahakis fc Koniglll200l . The same decrease in electron 
number that causes a drop in scattered X-ray emission 
near the disk would lead to a time-delayed increase in 
the speed of the jet downstream. The increase in flow 
speed of the jet could form a sho ck wave, eventuall y 
seen as a superluminal radio knot ([Gomez eFall [19971) . 
This would then give rise to the dip-ejection sequence 
discussed above. 

The highest-amplitude 37 GHz outburst, from 2007.5 
to 2008.8, started at the same time as at 230 GHz, but 
reached maximum level at a later time than both the 
230 GHz and X-ray flares. This can be explained by 
the larger optical depth at the lower frequencies prior to 
the 37 GHz peak. The similar amplitudes of the 37 and 
230 GHz flares follow the patter n of the synchrotron-los s 
stage of the shock- in-jet model (iMarscher fc Gea"i1ll985l ) , 
which predicts that, in synchrotron flares produced in 
the jet, the peak amplitude should stay roughly constant 
as the emission becomes optically thin at progressively 
lower frequencies. 

The mean time delay between the minimum of 
the X-ray dips and the time of ejection of the cor- 
responding superluminal knot is 0.15 ± 0.08 yr. The 
average apparent speed of the moving components 
with well-determined motions is 3.9c ± 0.7c. Therefore, 
a knot moves a distance of 0.18 ± 0.1 pc in 0.15 yr, 
projected on the plane of the sky. (Here we assume that 
acceleration of the flow to its terminal velocity occurs 
over a sufficiently short distance that the time that a 
moving feature spends in the acceleration zone is short 
compared with the total transit time to the 43 GHz 
core.) Since the angle of the jet axis of 3C 111 to the 
line of sight ~18°, the actual distance traveled by the 
knot, given by l3c6Atohs^ is 0.6 ± 0.3 pc. Hence, we 
derive a distance 0.6 ± 0.3 pc from the corona or the 
base of the jet (where the X-rays are produced) to the 
core seen on the 43 GHz VLBA image s. This is one of 
the fe w cases where, similar to 3C 120 ([Chatteriee et al.l 
l2009h . we are able to specify the distance between the 
central engine and mm-wave core in an AGN. 
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Figure 17. Discrete cross-correlation function (DCCF) of the 
optical and X-ray monitor data of 3C 111 for the entire 5 yr 
interval. The time delay is defined as positive if the X-ray 
variations lag those at the optical frequency. The jagged solid 
and dotted lines at the top of the panel denote 99 and 95 percent 
significance levels of correlation (see text), respectively. Those at 
the bottom of the panel show the same for anti-correlation. The 
peak X-ray versus optical DCCF is 0.6, which corresponds to more 
than 99% significance level. The peak at —15 days indicates that 
the X-ray variations lead those in the optical by 15 ± 10 days. 

6. X-RAY/OPTICAL CORRELATION 

As determined by the DCCF (Figure fT7|) . we find that 
the X-ray variations of 3C 111 are very strongly corre- 
lated with those at optical R band. The peak X-ray 
versus optical DCCF is 0.6, which corresponds to more 
than 99% significance level. The jagged solid and dotted 
lines at the top of the panel denote 99 and 95 percent sig- 
nificance levels of correlation, respectively. Those at the 
bottom of the panel show the same for anti-correlation, 
as described in the previous section. The time lag of the 
peak indicates that the X-ray variations lead those in 
the optical by 15 ± 10 days. The highly significant corre- 
lation between the X-ray and optical variations implies 
that emission at these wave bands is causally connected. 

Fig. [S] shows the variation of the degree of polarization 
as well as the electric vector position angle (EVPA) of 
3C 111 at R band from 2006 to the present. The av- 
erage polarization was only 1.6% ± 0.6%, rather low for 
synchrotron radiation from the jet. This agrees with the 
idea that most of the optical emission is th ermal radia- 
tion from the accretion disk (|Malkanlll983l ). The strong 
X-ray/optical correlation, weak optical polarization, and 
smaller variance of optical compared with X-ray flux 
at shorter timescales are consistent with a reprocessing 
model. In this scenario, the X-rays are predominantly 
produced by inverse Compton (IC) scattering of ther- 
mal optical/UV seed photons from the accretion disk by 
hot, but non-relativistic, electrons in the corona, while a 
significant fraction of the optical-UV emission is due to 
heating of the accretion disk by X-rays produced in the 
above process. The slope of the optical PSD in the range 
lO"^'"'^ to 10~^'* Hz is '^—2.0, significantly steeper than 
that of the X-ray PSD over the same frequency range (~ 
— 1.0). This implies that the optical variations on shorter 
time-scales are suppressed, consistent with smoothing by 
reprocessing. 

Alternatively, a substantial fraction of the optical emis- 



sion may be generated as synchrotron radiation in the jet, 
giving rise to a degree of polariza tion up to 3% (Fig. [SJ or 
even higher ([Jorstad et al.|[2007l) and significant variabil- 
ity of the position angle of polarization (between 0° and 
150°) at some epochs. In this case, if X-rays are indeed 
generated at the base of the jet, then the highly signif- 
icant correlation between the X-ray and optical varia- 
tions, combined with the time delay of the latter, 15 ± 10 
days, implies that the optical emission region in the jet, 
if present, is situated downstream of the base of the jet. 
We can calculate the mean distance (Az)xo between the 
X-ray and optical synchrotron emission sites in the jet 
if we adopt 18° as the angle betwe en the jet axis and 
the line of sight ([Jorstad et al.|[2005l ) and assume a value 
for the average speed f3c of the jet flow between the two 
emission regions: (Az)xo ~ 6, 20, and 200 light-days for 
/3 = 0.25, 0.5, and 0.974, respectively. The last value, 
/3 — 0.974, corresponds to the velocity obtained from 
analysis of the superluminal apparent motion and time- 
scale of flux decli ne at 43 GHz of knots observed with 
the VLB A (Jorsta d et all 120051) . The lower values of 13 
near the base of the jet are in concert with theoretical 
models in which jets are accelerated to relativistic flow 
veloc ities over extended distan ces from the central engine 
(e.g., IVlahakis fc K6nigll[200l . 

Some of the thermal optical-UV emission could be pro- 
duced in the accretion disk even without reprocessing 
via X-ray heating. The temporal nature of this emis- 
sion would be different from that generated as a result 
of reprocessing. The X-rays may be produced mainly by 
up-scattering of UV rather than optical photons. This 
could occur if the flux of optical photons reaching the 
corona is much smaller than that of UV photons. Such 
a scenario is likely if the corona is small, such that the 
UV emission region is much closer to the corona than 
the region where t he majority of the op tical photons 
are produced (see iChatteriee et al.l 120091 for details). 
Any disturbance propagating outward in the accretion 
disk will cause a change in the UV flux (and a resul- 
tant nearly immediate change in the X-ray emission) 
followed by a similar change in the optical fiux. This 
may also give rise to the X-ray/optica l time delay that 
is obs erved. As discussed in detail by IChatteriee et al.l 
(|2009| ) , the variability timescale describ ed above i s con - 
sist ent with the mod el pr oposed by K ing et al.l ()2004f ) 
and lLivio. Pringle fc Kind (j2003D . in which variability in 
the disk emission is caused by large-scale alignment of 
poloidal magnetic field in the inner accretion disk from 
random fluctuations in field direction. Such alignment 
occurs on a timescale of a few tens of days for a BH of 
mass of '^IO^Mq, consistent with the results described 
above. 

The time delay of 15 ± 10 days of the optical with 
respect to the X-ray variations is larger than expected 
from a pure reproces sing model. As discussed in 
iChatterjee et al.l ()2009f ). the region of the accretion disk 
that produces the largest amount of direct (rather than 
reprocessed) optical emission is at '^100 gravitational 
radii (r^) from the center, which is equal to ^1 light- 
day for a B H mass of 1.8 x 10^ M(t, (see a lso the detailed 
modeling of iKaz anas fc NayakshinI 120011 ). Therefore, if 
the X-ray/optical time delay is solely due to light travel 
time from the corona to the region in the accretion disk 
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producing reprocessed optical emission, the correspond- 
ing reprocessing region has to be farther and/or extended 
than this by an order of magnitude (^1000 r^). If, on 
the other hand, the variabihty is caused by disturbances 
propagating outward through the disk at velocity Vdist, 
then the optical emission region could be closer by a fac- 
tor Vdist/c. If a fraction of the optical emission is pro- 
duced in the jet, as discussed above, or the variability is 
partly affected by disturbances propagating in the accre- 
tion disk, then the observed time delay and flat peak of 
the X-ray/optical cross-correlation function could result. 

The accretion rate estimated in §3, m = 
(0.02 ± 0.01)modd, is within the range at which 
the transition from high to low sta te occurs in the 
Galactic BHXRBs (jMaccarond I2003D . An advection 
dominated accretion flow (ADAF) soluti on may exist fo r 
accretion onto black holes at this rate (|Naravanll2005D . 
The radiation in an ADAF is dominated by thermal 
comptonization, which has been used to explain the 
X-ray spect ra of BHXRBs in the canonical low-hard 
state (e.g., lEsin et all I1998D . It is possible that the 
X-rays from 3C 111 that we have attributed to a 
distribution of hot electrons consistent with the corona 
might instead be produced by an ADAF. 

7. COMPARISON WITH THE FR I RADIO GALAXY 3C 120 

Although 3C 111 and 3C 120 fall into two different 
Fanaroff-Riley categories, they have a number of com- 
mon characteristics. Both are radio-loud AGNs with jets 
that dominate the emission at radio frequencies, but with 
less radiative power than typical blazars and directed at 
a wider angle to our line of sight. In both radio galax- 
ies, the low optical polarization is consistent with most 
of the optical continuum resulting from blackbody radi- 
ation from the accretion disk. The X-ray properties are 
consistent with inverse Compton scattering of thermal 
optical-UV photons from the disk by a distribution of 
hot electrons — the corona, which could be the base of 
the jet — situated near the disk. In 3C 111, some of 
the optical and X-ray emission could be from the jet. A 
major difference between these two radio galaxies is in 
the accretion rate. In 3C 120, the accretion rate is high 
(~0.3 riicdd), similar to Galactic BHXRBs in the X-ray 
high state. In 3C 111, it is lower by an order of magnitude 
('-~^0.02 modd) which is close to the transition from high 
to low state in galactic BHXRBs. Therefore, the inner 
accretion flow in 3C 111 may be advection dominated. 

The X-ray emission from each of these radio galaxies 
possesses a characteristic timescale that is proportional 
to the BH mass, similar to what is observed in Galactic 
BHXRBs and Seyfert galaxies. This result implies that 
there is a universality in the accretion processes of BHs 
in the mass range IOMq to 1O^M0. These objects also 
show a clear connection between the radiative state near 
the BH and events in the jet. Ejection of new knots in 
the jet following dips in the X-ray light curve observed 
in these two radio galaxies is somewhat similar to the 
Galactic micro-quasar GRS 1915-1-105, where knots with 
apparent superluminal motion appear in the jet after the 
X-r ay flux recovers from intervals of very low flux lev- 
els ([Mirabel et al.l 119981 ) . This connection is expected 
according to current theories of jet launching and col- 
limation, but has not been demonstrated previously to 



extend to radio-loud AGNs. 

The relationships that we have uncovered between 
the central engine and the jet p rovide strong support 
for the paradigm (e . g.. iKprding. Jes ter & Fender 200(| 
iMarecki fc Swobodal l2011h that AGNs and Galactic 
BHXRBs are fundamentally similar, with characteristic 
time and size scales proportional to the mass of the 
central BH. This implies that we can develop and 
test models of AGNs based partly on observations of 
BHXRBs. 

8. SUMMARY AND CONCLUSIONS 

This paper presents well-sampled, 6-yr-long light 
curves of 3C 111 between 2004 and 2010 at X-ray, optical, 
and radio wavebands, as well as monthly images obtained 
with the VLBA at 43 GHz. The appearance of new knots 
in the VLBA images of the jet follow prominent dips in 
the smoothed X-ray light curve, demonstrating a con- 
nection between events in the jet at parsec scales and 
those occurring near or in the accretion disk. We have 
calculated the broad-band X-ray power spectral density 
(PSD) of 3C 111 and cross-correlated the X-ray and op- 
tical light curves, determining the significance of the cor- 
relations with light curves simulated in accordance with 
the PSD. From the analysis of the temporal properties of 
the observed variability in multiple wave bands, we infer 
the locations and radiative mechanisms of the X-ray and 
optical emission. 

Our main conclusions are as follows: 

(1) The X-ray PSD of 3C 111 is best fit by a bending 
power-law model with a smooth change in slope above 
a break frequency. The best-fit value of the break fre- 
quency is IQ~^ 05 jj2, corresponding to a break timescale 
of 13 days, which agrees with the empirical relation 
between bre ak timescale, BH mass , and accretion rate 
obtained by iMcHardy et al.l (|2006[ ) based on a sample 
of mainly radio-quiet Seyfert galaxies. This indicates 
that the accretion process in 3C 111 is similar to that 
of BHXRBs and Seyfert galaxies, and that there is a 
possible universality in the accretion processes in BHs 
spanning a broad range of masses from 10 Mq to 10^ 

Mq. 

(2) In 3C 111, the 2-10 keV continuum flux and the 
Fe line intensity are strongly correlated, with a time 
lag shorter than 90 days and consistent with zero. This 
implies that the Fe line is generated within 90 light-days 
of the source of the X-ray continuum. 

(3) All X-ray dips are followed by the ejection of a 
new knot in the VLBA images. We derive a distance 
from the base of the jet or corona (where the X-rays 
are produced) to the VLBA 43 GHz core region of 
~0.6 pc using the average time delay between the 
minimum of the X-ray dips and the time of passage of 
the corresponding superluminal knots through the core. 

(4) The dip-ejection connection can be explained if a 
decrease in the X-ray production is linked with increased 
speed in the jet flow, causing a shock front to eventually 
form and move downstream. This property of 3C 111 
is similar to 3C 120 and reminiscent of Galactic BH 
systems, where transitions to high-soft X-ray states are 
associated with the emergence of very bright features 
that proceed to propagate down the radio jet. 

(5) Most of the optical emission in 3C 111 is generated 
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in the accretion disk due to X-ray heating. However, 
weak optical polarization indicates that some fraction of 
the optical emission is synchrotron radiation from the jet. 
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